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ABSTRACT 

A detailed analysis of the nonthermal X-ray emission from the North- Western 
and Southern parts of the supernova remnant (SNR) HESS J1731— 347 with 
Suzaku is presented. The shell portions covered by the observations emit hard 
and line-less X-rays. The spectrum can be reproduced by a simple absorbed 
power-law model with a photon index F of 1.8-2.7 and an absorption column 
density Nn of (1.0-2.1) x 10^^ cm~^. These quantities change significantly from 
region to region; the North- Western part of the SNR has the hardest and most 
absorbed spectrum. The Western part of the X-ray shell has a smaller curvature 
than North- Western and Southern shell segments. A comparison of the X-ray 
morphology to the Very High Energy (VHE) gamma-ray and radio images was 
performed. The efficiency of electron acceleration and emission mechanism in 
each portion of the shell are discussed. Thermal X-ray emission from the SNR 
was searched for but could not be detected at a significant level. 
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Subject headings: acceleration of particles — ISM: individual (HESS J1731— 347) 
— X-rays: ISM 



1. Introduction 



Supernova remnants (SNRs) are widely believed to be the main source of Galactic 
cosmic rays. Observationally, several SNRs are found to have X-ray synchrotron-emitting 
shells, and are therefore widely known as electron accelerators at least up to ~TeV energies. 
Many of them are identified as young SNRs from the records of superno vae, expansion 
veloci ty , and their bri g ht ther mal X-rays (e.g. Gas A, Tycho, Kepler, SN 1006: lKoyama et al. 



(Il995l ): iBamba et al.l (l2005al )). However, the nature of the prominent non-thermal X-ray- 
emitting SNRs RX J 171 3— 3946 and V ela Jr. is still largely unclear (IKoyama et al.l 119971 : 
Slane et al.ll200ll : iBamba et al.l l2005bl ). The very low level of radio synchrotron emission, 
the lack of thermal X-ray emission, and poor expansion velocity data have yielded only 
weakly constrained key parameters of the SNRs, such as age and nature of the progenitor 
star. Remarkably, these SNRs are bright Very High Energy (V HE) gamma-ray emitters 
compared with other young SNRs (lAharonian et al.l l2004j . |2005| ). Thus, they are strong 
particle accelerators, and it is important to understand the nature of this small class of 
sources, and especially the puzzling absence of thermal X-rays. 

HESS J1731— 347 has b een discovered in the VHE Galactic plane survey performed 
with the H.E.S.S. telescopes (lAharonian et al.ll2008l ). iTian et al.l ( 120081 ) discovered a shell 
structure in archival radio data spatially coincident with the VHE source, and the source 
(also named G353.6— 0.7) was hence considered to be a newly discovered SNR with associated 
VHE gamma-ray emissio n. Subsequently, deep VHE follow-up observations were performed 
(lAbramowski et al.ll201lt ) that established a clear VHE shell-like structure in agreement with 
the radio shell. 



Based upon a possible association of the SNR with a nearby HII region, iTian et al. 



(120081 ) suggested that the object could have the same distance (3.2 ± 0.8 kpc), using a 
fiat Galactic rotation c urve model and the most recent estimates of the parameters by 
Eisenhauer et al.l (120051 ). A lower limit of the distance comes from the clear correlation 
of the X-ray abs orption column density w ith a molecular cloud structure at ~3 kpc as seen 
in GO emission (lAbramowski et al.ll201ll ). suggesting that the cloud is in the foreground of 
or interacting w ith the SNR. This SNR has no GeV counterpart in the Fermi 2-year Point 
Source Gatalog (IFermi-LAT collaborationll201ll ). 
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XMM-Newton and Suzaku observations towards the center and the East of this SNR 



show that X-ray emissi on from the shell is dominated by non-thermal emission ( lAcero et al. 



20091 : iTian et al.ll2010f ). This suggests that HESS J1731-347 could be the third member 



of the SNR class characterized by pure non-thermal X-ray emission and comparable VHE 
gamma-ray luminosity. A sensitive search for thermal X-ray emission across the entire rem- 
nant is needed to confirm this hypothesis. 

From spatially resolved X-ray synchrotron spectra, changes of particle acceleration con- 
ditions can be studied. Such variations have already been found in other SNRs such as 
RX J1713 — 3946, and were att ributed to electron cooling over the distance from the acceler- 
ation site (ITanaka et al.ll2008l ). 



In this paper, we report on Suzaku observations (IMitsuda et al.l 120071 ) of the western 
half of the SNR HESS J1731— 347. The observation details are summarized in §2J A first 
analysis of the Suzaku X-ray data is presented in ^ the results of which are discussed in §11 



Observations and Data Reduction 



Suzaku observed the center region of HESS J1731— 347 on 2007 Feb. 23-24, and the 
North- Western and Southern regions on 2010 Feb. 17-19. The mapping consists of three 
pointings as listed in Table [1] Data reduction and analysis were made using HEADAS 
software version 6.9, version 2.5.16.28 of the processed data, and XSPEC version 12.6.0. 

Suzaku has tw o active instruments, four X-ray Imaging Spe ctrometers (XIS0-XIS3 ; 
Koyama et al.ll2007l ) each at the focus of an X-Ray Telescope (XRT s ISerlemitsos et al.ll2007l ) 
and a separate Hard X-ray Detector (HXD; iTakahashi et al.l 120071 ). 



Only three XISs could be used subsequently due to a problem with XIS 2. XISl is 
a back-illuminated (BI) CCD, whereas the others are front-illuminated (FI). The XIS in- 
struments were operated in normal full-frarn e clocking mode. Spaced-row charge injection 
(jNakajima et al.l l2008l : lUchiyama et al.l 120091 ) was used in the later two observations (Ta- 
ble [1]). We filtered out data acquired during passages through the South Atlantic Anomaly 
(SAA), with elevation angle to the Earth's dark limb below 5°, or with elevation angle to 
the bright limb below 25° in order to avoid contamination by emission from the bright limb. 
The remaining exposure time for each observation is summarized in Table [TJ 

The HXD PIN was operated in normal mode. We filtered out data obtained during 
passages through the SAA, with elevation angle to the Earth's limb below 5°, and cut off 
rigidity (COR) smaller than 8 GV. The exposure time for each observation is again listed 
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in Table [H For the no n- X-ray background (NXB) model, we adopted the LCFIT model by 



Fukazawa et al. (12009 ) . The cosmic X-ray background (CXB) flux is estimated based on the 



HEAOl results (iBoldtl 119871 ). and treated as an additional background component. 



Results 



3.1. X-ray images 



Figured] shows 0.5-2.0 keV and 2.0-8.0 keV mosaic images obtained with the Suzaku XIS 
instruments. We used only XISl and XIS3 for image analysis, since a part of XISO showed a 
known anomaly (JX-ISAS-SUZAKU-MEMO-2010-01) and discriminated that region in the 
data pro cess before deriving spectra. The vignetting has been corrected in each image using 
xissim (jlshisaki et al.l 120071 ). 



One can see clear shell-like structures in both images. The shell has a rather clumpy 
and bulged structure. The Western part of the X-ray shell has a smaller curvature than 
the North- Western and Southern shell segments. At the center of the remnant, emission 
from the point source XMMU J173203. 3— 344518 (ITian et al.ll2008l ) is seen, which from its 
soft spectrum and its position has been c l assified as the cen tral compact object (CCO) 
(lAcero et al.ll2009l : lAbramowski et al.l 1201 ll ). iTian et al.l (|2010( ) noted that the source spec- 
trum is also compatible with a magnetar hypot hesis. A claim of X- r ay pul sation with a 
period of ~1 s using XMM-Newton data also led iHalpern &: Gotthelll (l2010a| ) to speculate 
that XMMU J173203. 3— 344518 might be a magnetar candidate, but a subsequent dedi- 
cated search for pulsations using Chandra timing observations did not confirm the signal 
faalpern fc Gottheljl2010bl ). 



In Fig. 121 a false-color image of the remnant is shown. The image reveals that the 
North- Western part of this object is characterized by significantly harder emission compared 
with the other regions. 



3.2. Multiwavelength comparison 



This SN R has a shell-like structure not only in X-rays but also in radio and VHE 
gamma-rays ( iTian et al.l l2008l : lAbramowski et al.l 1201 ll ) . The top panel of Figure E] shows a 
comparison of the Suzaku data with a ra dio map at 1.4 GHz from the ATCA south Galac- 
tic Plane survey (iHaverkorn et al.ll2006l ). Both images show clumpy shell structures, and 
an overall good agreement with each other. The bottom panel of Fig. [3] shows a compari- 
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son of the Suzaku data y ith the VHE gamma-ray map derived from H.E.S.S. observations 
( lAbramowski et al.ll201l[ ). At first glance, it looks as if there are differences in the extent of 
the VHE and X-ray emissions towards the North- West and West of the source, respectively. 
Those differences might however be attributed to the coarser angular resolution of the VHE 
data (0.06 deg., 68% con tainment radius) com pared with the Suzaku X-r ay data (half power 



diame ter of ~2 arcmin, iMitsuda et al.l 120071 ) . We refer to the result of lAbramowski et al. 



(120111 ) that the radio and VHE radio profiles of the entire remnant match well when the 
radio data are blurred to VHE resolution. 

To illustrate the potential difference between the X-ray and the VHE radial profiles, 
in Fig. m the 2-8 keV profile for a wedge towards the west is compared to the VHE profile 
derived from the entire remnant. Besides the difference in angular resolution, it is also 
not yet feasible with the current VHE statistics to restrict the VHE data to a similarly 
small region. Further V HE observations using e.g. the planned Cherenkov Telescope Array 
CTA (lActis et al.l 1201 ll ) are needed to provide sufficient VHE photon statistics to permit 
morphological comparisons on such small angular scales. 



3.3. Spectra below 10 keV 

3.3.1. Diffuse emission 

In order to perform a spatially-resolved spectral study below 10 keV with XIS, we 
divided the X-ray map into several regions, as shown in Figure [H The background regions 
were taken from source-free regions of the same observations. 

As can be seen in Figure [5l all background-subtracted spectra are hard and show no 
evidence of line emission. We thus fitted the sp ectra with an absorbed power-law r nodel . 
The absorption model i ncludes the cross sections of lBalucinska-Church fc McCammonI (119921 ) 
with solar abundance (lAnders fc Grevessd Il989l). Auxiliary files for the effective area were 



produced with xissimarfgen (llshisaki et al.ll2007l ). using the observed images of the source 
regions as photon distribution data. For each region, the spectra of all available XIS detectors 
were fitted simultaneously. For some regions, less than three XIS spectra are available 
because the regions overlap with the damaged part of XISO or with calibration sources. 
Because of this, we allowed for a difference of the normalization between different XIS spectra, 
and adopted the best-fit normalization derived from XIS3 as final result, since any remaining 
contamination from calibration sources is smallest in XIS3 for the regions analysed. Figure E] 
and Table [2] show our best-fit models and parameters with 90% errors. In each case the 
power-law fit was statistically acceptable. 
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Figure [8] shows the best-fit photon index (left) and column density (right) for all analysed 
regions of the SNR. The North- Western part of the shell has a significantly harder and 
more absorbed spectrum, compared with the rest of the covered remnant. On the other 
hand, differences between other individual regions are not significant when statistical and 
systematic errors (at low surface-brightness regions, see below) are taken into account. 

In order to better constrain the spectral shape of the non-thermal emission, we combined 
several adjacent regions to three cardinal areas: regions 1, 4, and 5 to the North- West, regions 
3 and 6 to the Center-East, and regions 8-13 to the South, respectively. All the three spectra 
were again fitted well with an absorbed power-law model without any spectral break. The 
spectra with the best-fit models and their best-fit parameters are shown in Fig. |5]and Tab. [2l 
respectively. As already seen from the individual spectra, also in the combined spectra there 
is a significant change in photon index across the remnant, with the hardest spectrum in 
the North- West, where also the absorption column density is significantly larger than in the 
other two 

As illustrated in Fig. [H we used two apparently source-free regions to estimate the aver- 
age background of the entire source region. Due to the position dependence of t he Galactic 



Ridg e X-Ray Emission (GRXE), that can be ~20% contribution in this region (jUchiyama 



20101), the spectral parameters of low surface brightness regions might be misestimated. 
In order to quantify this systematic error, we artificially varied the background by ±20% 
in the fit of region 9, which has the lowest surface brightness of all regions analyzed. In- 
deed, the systematic error in this region is comparable with the statistical error: F = 3.2 
(+0.5-0.5)stat(+0.7-0.4),ys and A^H = (2.5 (+0.7-0.6)stat(+0.9-0.5)sys) xlO^^ cm^. For the 
high surface brightness regions (> 0.8 x 10^^^ erg cm^^s^^arcmin^^, see Table |2]) which are 
relevant for the evaluation of the spectral and absorption trends, the systematic errors are 
below the 90% statistical errors. Therefore, while the significance of the observed absorption 
and photon index trends indeed weakens, they are nevertheless significantly observed. 

In a power-law spectrum, photon index and absorption column densities are correlated 
to some extent. To further verify that the change of F and A^h in the Suzaku data across 
the remnant is indeed significant, we compared the confidence contours of F and A^h in the 
three cardinal areas. The left panel of Fig. [7] confirms that the North- Western area has 
significantly harder emission and higher absorption than the other two areas, and the south 
area has the softest emission. Similar checks were also done for several individual regions; 
the right panel of Fig. [7] illustrates that region 5 indeed exhibits significantly harder and 
more absorbed emission than region 12. 



To investigate whether the source emits thermal radiation, we searched for emission 
lines of thermally heated plasma on top of the continuum spectra, but could not identify 
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promising candidates. This applies to the spectra derived from individual regions, as well 
as to the combined spectra from the three cardinal areas discussed above. We also tested 
the possibility that the lineless spectra could come from a highly absorbed hot plasma with 
low ionization age instead of from a nonthermal electron plasma. For this purpose, we fitted 
the softest spectrum from region 9 which also has the highest absorption column density, 



with an absorbed non-equilibrium ionization plasma model (NEI model, iBorkowski et al. 



2OOII ). Indeed, the fit is statistically acceptable with a reduced of 47.2/39. The resulting 



temperature is however higher than ~2 keV, with an ionization timescale of ~ 10^ s cm~^. 
Same conclusion or result holds if several regions with soft spectra are combined, using e.g. 
regions with photon indices larger than 2.5 in the sourthern half, namely regions 7, 8, 9, 
11, 12, and 13. Again, the NEI model reproduces the spectra well with a reduced of 
60.7/118, with similar parameters as in region 9 [kT >2.4 keV, ionization time scale less 
than ~ 2 X 10^ s cm~^). Such parameters would require a very young SNR or a recent 
encounter of a fast shock with a dense region, e.g. in a wind bubble scenario. Since age 
and environment of HESS J1731— 347 are not known, such an interpretation is formally 
not precluded. It is however more natural to assume that the spectra in these regions are 
dominated by nonthermal synchrotron emission, in agreement with the hard regions where 
a pure thermal interpretation is virtually excluded. 



3. 3. 2. Spectrum of XMMU Jl 73203. 3- 344518 

The central source of the SNR, XMMU J173203.3-344518, is in the field of view of 
the Suzaku XISs only in the first Suzaku pointing (ObsID 401099010, 2007/ 02/23-24) 



A firs t analysis of that data from the object has already been published by iTian et al. 



(120101 ). The uncertainty of their analysis is from the combination of the relatively coarse 
angular resolution (~2 arcmin) of Suzaku and the contamination from the surrounding SNR 
emission. We therefore repeated the spectral extraction, adding systematic checks to verify 
the robustness of the spectral results. 

We used a circular source extraction region of 2.5 arcmin radius centere d on XMMU J1732 03.3— 34451 



for th e spectral analysis, which encompasses ~90% of all source photons (ISerlemitsos et al. 



20071 ). As nominal background region, we chose an annulus around the source with 3-4.5 ar- 



cmin radius. 

As already known, the central source has a rather soft spectrum compared with the 
diffuse emission from the SNR. As the first assumption, an absorbed blackbody model and 
an absorbed power-law model were fit to the data. The fits are statistically unacceptable 
(reduced of 388.8/304 or 443.8/304, respectively). We conclude that the spectrum has a 
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multi component emission. 

Fitting with a blackbody plus a power-law component model or a two-temperature 
blackbody model, yields better results, with reduced of 334.7/302 or 340.7/302. The 
bes t-fit param e ters fo r the blackbody plus power-law model are consistent with the results 



by iTian et al.l ( 20101 ). The photon index for the power-law component, 4.7!!^q5, rather 
soft for magnetars and CCOs; we therefore used a two temperature blackbody model in this 
paper. We show the spectrum with the best-fit model in Fig. E] and the best-fit parameters 
in Tab. [31 

The second component, on the other hand, could be the contamination of the SNR 
emission. To check this possibility, we also fitted the spectrum with absorbed power-law or 
blackbody model plus the spectrum of the region 3 as a contamination, which is in the vicinity 
of the central source. The normalization of the contamination component was treated to be 
a free parameter, whereas photon index and absorption column density are fixed to those for 
region 3 spectrum. The x^/d.o.f. becomes 443.8/303 for the power-law case and 350.8/303 
for the blackbody case. The latter case has similar x^/d.o.f. as the multi-component models. 
It is found that the source spectrum can be well reproduced by a single blackbody model 
and the second component in the spectra could be just the contamination of the diffuse SNR 
emission. The best-fit parameters of the previous model are shown in Tab. |3l 



3.4. Spectra above 10 keV 



We also analysed the HXD PIN data-set to search for a signal from the source in 
the energy range above 10 keV. After background (NXB-I-CXB) subtraction, the remain- 
ing count rate in the 15-40 keV band in each observation is 3.7 ± 0.3 x 10^^ cuts s~^ for 
OBSID=401099010, 1.7 ±0.3 x 10"^ cuts s-^ for OBSID=504031010, and 2.3 ± 0.3 x lO'^ 
cuts s-i for OBSID=504032010. This is 12.4%, 2.9%, and 7.9% of the NXB count rate in 



this band, respectively. The systematic NXB uncertainty is 3-5% (IFukazawa et al.l 120091 ). 
i.e. ~1-I.5xl0~^ cuts s~^, b ut the uncertainty cou ld be larger when the observation du- 
ration is shorter than 1 day (IFukazawa et al.l |2009| ) as is the case here. Ignoring for the 
moment this caveat, the first and third observation appear to display significant emission 
above instrumental and extragalactic backgrounds in the 15-40 keV background. 

Before potentially being able to attribute the excess signal to HESS J1731— 347, another 
background, the GRXE, has to be considered. Its level is difficult to estimate sinc e it is rather 
faint a bove 10 keV and its measurement is contaminated by bright point sources. iRevnivtsev 
( 12003! ) revealed with RXTEXh-oX the GRXE in this region is a few % of the diffuse emission in 
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the Galactic center (GC) region. This result is consistent with results obtained with GINGA 
X-ray satellite (Yamauchi and Ebisawa, private communication). With HXD PIN, the GC 
diffuse emission was resolved to be 4 x 10~^° erg cm~^s ~^deg:~^ in the 12-40 keV band, or 
~ 0.5 cnt s~^ in the 10-40 IceV band (lYuasa et al.ll2008l ). Our count rate is roughly 10% of 
that in the GC region, thus the emission is marginally above the GRXE. 

Since HXD has no spatial resolution, any point source in the field of view (~0.5 deg. 
FWHM for PIN) adds to the total count rate in the detector. First of all, the CCO, 
XMMU J173203.3-344518, is almost at the center of OBSID=401099010. However, its 
spectrum is very soft, well described by a thermal blackbo dy with temperature 0.5 keV 
(lAcero et al.l 120091 : iTian et al.ll2010l : lAbramowski et al.ll201l[ ). and does therefore not signif- 
icantly contribute to the emission above 10 keV in the field. 

To identify potential hard-band emitting sources, we compared the three Suz aku fields of 



view with the INTEGRAL Fourth IBIS/ISGRI Soft Gamma-ray Survey Catalog (IBird et al. 
20101 ). As a result, we found 4U 1728—337, a low-mass X-ray binary, which is ~50 arcmin 



away from the center o f the SNR. The so urce is variable with a 40-70 keV flux of ~(1- 
6)xl0^^° ergs s~^cm^^ (jClaret et al.lll994l ). During our Suzaku observations, 4U 1728—337 
was detected with the RXTE All Sky Monitor (ASM) with a count rate of a few cuts s^^. 
Converting this value into an on-axis HXD PIN count rate results in a rough count rate 
estimate of ~0.1-1 cuts s~^ in the 15-40 keV band. For the actual source location at the 
edge of the HXD PIN FOV the expected count rate is one order of magnitude smaller, 
i.e. ~0. 01-0.1 cuts s^^. Given that this count rate is at the same level as the measured 
excess, from our observations we cannot claim any detection of X-r ays above 10 keV frorn 
HESS J1731-3 47. Further studies using forthcoming ASTRO-H jTakahashi et all boioh 
and/or NuStar (IHarrison et al.ll2010l ) are needed to confirm this. 



3.5. Searching for coherent pulsations from the central compact object 



X-ray pulsations in the range of 0.1-0.4 s have been detected from a few CCOs (iHalpern fc Gotthelf 



2010d ). On the other hand, magnetars typically ex hibit pulsations in the rang e of 2-12 sec- 
onds. The fact that the claimed pulsa tion of ~1 s (IHalpern fc Gotthelfll2010al ) could not be 
confirmed jHalpern fc Gotthelj[2010b[ ) does not preclude pulsations from XMMU J173203.3-344518. 



Below 10 keV, the time resolution of XIS is 8 s (IKoyama et al.l 120071 ). insufficient for 
a sensitive search for pulsations. Above 10 keV, HXD provides adequate time resolution 
(61 /is ), and systenaatic ba ckground effects are negligible when searching for coherent pulsa- 
tions JTerada et al.ll2008al lbll Hence, even though the source was not significantly detected 
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in the HXD spectral analysis, a pulsation search was still performed. However, no significant 
pulsations were found in the scanned period range of 0.1 s to 1000 s. 



Discussion 



4.1. Diffuse emission from SNRs 



The detailed analysis of the Suzaku data revealed that the North- Western shell of 
HESS J1731— 347 exhibits stronger absorption than the rest of the source. The same 
trend has already been observed with XMM-Newton (lAbramowski et al.ll201ll ). Actually, 
the statistics obtained with Suzaku permit to test smaller scales, but the two regions (4 
and 9, see Fig. [H]) that seem to deviate from the overall trend also have the lowest surface 
brightness, and further investigation is needed to verify the result. 

Towards the North- West of the SNR lies the Galactic Plane. However, the increase of 
the X-ray absorption column density can most likely not be attributed to the average gas 
density increase towards the Galactic Plane, since the measured iVn (1.9-2.1x10^^ cm~^) 
is in clear exces s of the average expected absorption column density in this direction (1.2- 
1.3x10^^ cm~^ (jPickev Sz Lockmanlll990l : iKalberla et al.ll2005l )). That value is in fact consis- 
tent with the overa l l colu mn density measured towards the South-East of the source. Indeed, 
Abramowski et al.l (120111 ) showed that the absorption column derived from the XMM-Newton 
Eastern coverage of the source is well correlated with molecular gas, the density of which is 
derived from CO observations towards that region. The question remains whether the SNR 
is located at a distance of ~3 kpc and interacting with a dense molecular cloud identified 
through the CO measurement at that distance, or whether the source is in fact located be- 
hind the cloud. CO observations of higher excited states may be helpful to find shock-heated 
gas. 

The non- uniformity of the measured X-ray photon indices might offer another clue to 
that question. In SNRs, hard synchrotron spectra such as the ones detected here are typically 
attributed to electrons acceleration in fast shocks (several thousand km s~^). The harder 
the spectrum, the higher the maximum particle energy and hence the higher the speed of 
the shock. A lower shock speed might indicate the presence of dense material into which the 
shock is propagating. However, in the observations presented here, the North- West, which 
displays the highest column densities, also shows the hardest photon indices (e.g., region 1 
and 5 in Table [2]), at odds with a cloud-interaction scenario. We therefore stress again that 
at present there is no observational implication of the SNR-molecular cloud interaction in 
the North- Western part, and that it is more likely that the absorbing material lies in front of 
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the source. Furthermore, the flat azimuthal VHE profile and the lack of correlation between 
the absorption column and X-ray surface brightness also indicate that the SNR may in fact 
not be interacting with the molecular cloud but rather that the source is located behind the 
cloud. 

The accumulated unabsorbed X-ray flux in our analysed regions is 2.3 x 10~^^ ergs cm~^s~^ 
in the 2-10 keV band. For the entire remnant, usin g our Suzaku observations together with 
the XMM-Newton results ( lAbramowski et al.ll201ll ). the total flux of the source in this band 
is approximately ~ 5 x 10~^^ erg cm~^s~^. The non-thermal X-ray luminosity L and the 
radius of the source R are estimated to L ~ 5 x 10^^ erg s~^ x d^kpc' and i? ~ 14 pc x ciskpc 
(c?3kpc = d/{?) kpc), where d is the distance to the remnant). For a distance of 3 kpc, these 
values ar e in rough agree n aent w ith the SNR radius vs. non-thermal X-ray luminosity re- 
lation of iNakamura et al.l (120111 ). If this relation also holds for HESS J1731— 347, then it 
is possible that the SNR is in fact very near (but still behind) the molecular cloud that is 
responsible for the strong X-ray absorption. 

Typical photon indices in the soft X-ray band f rom young SNRs like RX J1713— 3946 
range between ~2.3 and ~2.6 ( iKoyama et al.l 119971 : ISlane et al.l 120011 ). In the North- West 
of HESS J1731— 347, F is as small a s ^1.8. The only SN R known so far which shows a 
comparable photon index is CTB 37B (INakamura et al.ll2009l ). though the spectral extraction 
in the region of interest was very complex and subject to systematic uncertainties. Recently, 
Takahashi et al.l (120081 ) discovered that the synchrotron spectra of RX J1713— 3946 have a 
spectral break in the hard X-ray band, and that the pho ton index below th e break is quite 
hard (~1.5). A similar break was also found in SN 1006 (IBamba et al.ll2008l ). We speculate 
that the hard spectra in the North- Western region of HESS J1731— 347 with F < 2 are of 
similar nature. The Western and Southern regions have softer X-r ay spectra. The pho ton 
index of 2.3-3.0 is typical for synchrotron X-rays from SNR shells (IBamba et al.ll2005al . for 
example). We note that the Western part of the X-ray shell has a smaller curvature than the 
North- Western and Southern shell segments. This area looks promising to evaluate whether 
the VHE emission here does follow this trend, or whether the VHE shell is compatible with 
the larger average shell radius as the images seem to indicate (see Fig. [3]). Such differences 
could be indicative of shock-cloud interaction and might provide important clues of shock- 
accelerated protons in SNRs: The X-rays trace the shock which has slowed down due to 
the encounter with a molecular cloud, while high energy protons have partially escaped 
the shock and are primarily emitting VHE 7-rays in th e dense molecula r cloud material. 
Similar phenomena are seen in W28 north-eastern shell (INakamura! |20 1 ll ) and G359. 1—0.5 
(IBamba et al.l 120091 ). As discussed earher in this paper, the available VHE statistics are 
presumably not sufficient for such a test at this point of time. Also, the coverage of the 
X-ray emission in this particular area should be extended towards the West to pin-point the 
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outermost X-ray boundary here. 

In several SNRs, synchrotron X-ray spe ctra are softer in down stream regions than 
regions close to the shock, e.g. in SN 1006 ( iRothenflug et al.l |2004| ) or in Tycho's SNR 
(jCassam-Chenai' et al.l 120071 ). Within the still limited Suzaku coverage of HESS J1731— 347, 
however, no systematic trend could be found in the photon indices when comparing regions 
which are apparently shock-dominated with regions closer to the center of the remnant. 



4.2. XMMU J173203.3-344518 



The spectrum of the central source XMMU J173203. 3— 344518 has at least one tem- 
perature blackbody component. If the hard component is represented by a power-law, the 
best-fit photon in dex, 4.7 (4.2-5.2), is not compatible with the hard tails usually detected 
in magnetars (c.f.. lEnoto et al.ll2010l ). Moreover, the best-fit absorption column of 2.5 (2.2- 
2.7) X 10^^ cm^ significantly exceeds the values measured from the diffuse emission of the 
SNR. We therefore adopt in the following the two temperature blackbody model for the 
emission from this source. 

Assuming that the distance to this source is 3 kpc, the radius of the emitting region is 
2.2 (1.4-4.4) km for the low kT component (= -Rlt) and 0.6 (0.2-0.8) km for the high kT 
component (= Rut), respectively, i.e. the radiation is emitted from localized spots on the 
neutron star surface. 

Although no flare has been detected from XMMU J173203. 3— 344518, it is instructive 
to compare the spectral results to flare spectra of soft gamma-ray repeaters (SGRs). These 
are well reproduced by two temperature blackbody models, a nd the radii of the two com- 
ponents has a positive correlation of /^ht^Z-Rlt^ ~0.03-0.4 (INakagawa et al.l 120071 ). For 
XMMU J173203. 3— 344518, the ratio is 0.06, consistent with the above relation. Therefore, 
although the best-fit temperatures of both components are lower than flares, it is suggestive 
to interpret the emission from XMMU J173203. 3— 344518 in a magnetar scenario. Never- 
theless, while a magnetar nature of XMMU J173203. 3— 344518 is possible, it clearly cannot 
be excluded, that the source is a CCO associated with the supernova remnant. Such objects 
are typically found in young SNRs and characterized by soft black body-like emission with 
kT ~ 0.5 keV and a lack of variability (except pulsations with periods below one second 
detected in some cases). The spectrum of one such object, PSR J1852-I-0040 — a confirmed 
CCO in SNR Kes 79, has also been fitted with a two -component black body {kT =0.3 and 
0.52 keV), as reported bv lHalpern fc Gotthelj J2010a[ ). 



On the other hand, the harder component can be just contamination of the diffuse 
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SNR emission. The source spectrum in this case, ~0.5 keV blackbody, is also typical for 
CCOs. Further detailed observations with good spatial resolution and statistics are needed 
to confirm the nature of this source. 
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Fig. 1. — Suzaku XIS intensity maps of the North- Western, Western, and Southern regions 
of HESS J1731-347 in the 0.5-2.0 keV (top) and 2.0-8.0 keV (bottom) bands. The image 
is binned with 4 arcsec and smoothed with a two dimentional Gaussian of cr = 12 arcsec. 
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Fig. 2. — A false-color reprensentation of the HESS J1731— 347 North- Western, Western, 
and Southern regions. Red and blue correspond to 0.5-2.0 keV and 2.0-8.0 keV X-rays, 
respectively. The image is binned with 4 arcsec and smoothed with a two-dimentional 
Gaussian of cr = 12 arcsec. 
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Fig. 3. — Multiwavelength view of the HESS J1731— 347 region. The color images show 
the Suzaku data in a false-color representation with the same color scaling as in Fig. Ej 
whereas the contours are derived from the A TCA radio map at 1.4 G Hz (top) and from the 
H.E.S.S. VHE gamma-ray excess count map (jAbramowski et ah 2011 ) (bottom). The green 
circle denotes the best-fit model of the VHE shell, which has a radius of 0.27 deg. and an 
unresolved shell width. 



- 19 - 



in 
d 



■ . HESS > 240 GeV, P.A. 0-360 deg 
o Suzaku 2-8 keV, P.A. 80-150 deg 



0.1 



0.2 0.3 
Radius (degrees) 



0.4 



Fig. 4. — VHE 7-ray (filled circle) and 2-8 keV (open circle) radial profiles. Both profiles 
are centered on the central compact object, (RA, Dec.) = (263.0125, -34.755). The central 
part was cut in the X-ray data due to the contamination from the central source. Note that 
the VHE 7-ray profile is taken from the entire remnant, whereas the X-ray one is from the 
Western part of the shell, using position angle (P-A.) 80°-150° (counted clockwize starting 
from North). 



Table 1. Observation Log 



ObsID 


Date 


Position 


XIS Exposure 


HXD Exposure 


SCI 




YYYY/MM/DD 


(J2000) 


[ksec] 


[ksec] 




401099010 


2007/02/23-24 


(263.0179, -34.7706) 


41 


35 


OFF 


504031010 


2010/02/18-19 


(263.0074, -34.9458) 


42 


34 


ON 


5010:-!-?()l() 


2010/02/17-18 


(202.8 18:-!, -:ii.()32."')) 


12 




ON 
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Fig. 5. — XIS spectra of each region (for the region definition see Fig. [T]) and combined one 
(see text). The color represents XISO (black), XISl (red), and XIS3 (green). The solid lines 
represent the best-fit absorbed power-law model. The lower panels show the residuals from 
the best-fit models, respectively. 
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Energy (keV) Energy (keV) 

Fig. 6.— XIS spectra of XMMU J173203.3-344518. The color represents XISO (black), 
XISl (red), and XIS3 (green). Left panel shows the two-temperature model fitting, whereas 
the right one shows a single blackbody plus SNR contamination model fitting. The lower 
panels show the residuals from the best-fit model. 




Fig. 7. — 68% (black), 90% (red), and 99% (green) confidence contours of the photon index 
vs. absorption column for North- West, Center-East, and South regions (left) and regions 5 
and 12 (right). 
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Fig. 8. — Best-fit photon index (left) and absorption column density (right) in each region. 



Table 2. Best-fit parameters of spectral fittings'^ . 



Region 








r 




F2 


- lOkcV 


Surface brightness 


xVd.o.f. 




(10 


cm~^) 








(10-12 


3rg cm^^s-i) 


(10" erg cm-2s-iarcmin-2) 




1 


2.11 


(1.92-2.30) 


2.21 


(2.09- 


-2.33) 


2.07 


(2.00-2.15) 


1.04 


(1.00-1.08) 


173.6/153 


2 


1.78 


(1.61-1.95) 


2.68 


(2.53- 


-2.84) 


1.66 


(1.61-1.95) 


1.38 


(1.34-1.63) 


98.5/93 


3 


1.40 


(1.34-1.47) 


2.37 


(2.31- 


-2.43) 


5.44 


(5.32-5.56) 


2.72 


(2.66-2.78) 


371.5/320 


4 


1.12 


(0.92-1.36) 


2.30 


(2.09- 


-2.54) 


0.73 


(0.67-0.78) 


0.42 


(0.38-0.45) 


111.7/79 


5 


1.93 


(1.73-2.15) 


1.77 


(1.64- 


-1.90) 


2.19 


(2.10-2.27) 


1.46 


(1.40-1.51) 


144.2/124 


6 


1.11 


(1.02-1.21) 


2.42 


(2.32- 


-2.53) 


2.01 


(1.94-2.09) 


1.68 


(1.62-1.74) 


137.8/136 


7 


0.96 


(0.82-1.12) 


2.59 


(2.38- 


-2.81) 


1.17 


(1.05-1.29) 


1.11 


(1.00-1.23) 


54.8/60 


8 


1.15 


(1.00-1.31) 


2.52 


(2.33- 


-2.71) 


1.54 


(1.43-1.66) 


0.88 


(0.82-0.95) 


119.0/94 


9 


2.52 


(1.96-3.18) 


3.17 


(2.70- 


-3.71) 


0.19 


(0.13-0.24) 


0.16 


(0.11-0.20) 


47.8/40 


10 


1.63 


(1.42-1.86) 


2.48 


(2.29- 


-2.68) 


1.56 


(1.43-1.69) 


1.27 


(1.17-1.38) 


49.7/74 


11 


1.10 


(1.00-1.20) 


2.54 


(2.43- 


-2.66) 


1.53 


(1.47-1.60) 


1.25 


(1.20-1.31) 


122.7/124 


12 


1.31 


(1.25-1.37) 


2.74 


(2.68- 


-2.81) 


3.20 


(3.12-3.27) 


1.31 


(1.27-1.33) 


380.4/337 


13 


1.40 


(1.18-1.65) 


2.55 


(2.34- 


-2.78) 


0.65 


(0.60-0.70) 


0.53 


(0.49-0.57) 


82.5/65 


North- West . 


1.66 


(1.54-1.79) 


1.97 


(1.88- 


-2.06) 


4.55 


(4.43-4.67) 


0.86 


(0.84-0.89) 


240.2/180 


Center-East 


1.29 


(1.24-1.35) 


2.36 


(2.30- 


-2.41) 


7.53 


(7.38-7.68) 


2.35 


(2.31-2.40) 


318.0/229 


South 


1.30 


(1.24-1.36) 


2.63 


(2.56- 


-2.70) 


7.74 


(7.50-7.98) 


0.85 


(0.82-0.88) 


347.4/355 




''Errors indicate single parameter 90% confidence regions. 



-23- 



Table 3. Best-fit parameters of spectral fittings of XMMU J173203.3-344518' 



Parameters 


2 blackbody case 


bla<;kbody + SNR contamination 


Nil (10^2 cm-2) 


1.7 (1.5-2.0) 


1.5 (1.4-1.6) 


fcTi (keV) 


0.33 (0.26-0.42) 


0.47 (0.46-0.48) 


Normi ^ 


14 (12-16) 




kT2 (keV) 


0.58 (0.53-0.72) 




Nor 1712^ 


3.5 (0.6-7.0) 




SNR contamination (10~^^ erg cm^^s^^arcmin"^)'^ 




5.6 (4.2-7.0) 



^Errors indicate single parameter 90% confidence regions. 

^ Rkm'^ / Dio^ , wiicre Rkm is the source radius in km, and, Dio is tlie distance to the source in units of 10 kpc. 
= 2-10 keV band. 



